Critical Josephson current through a bistable single-molecule junction 
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We compute the critical Josephson current through a single-molecule junction. As a model for 
a molecule with a bistable conformational degree of freedom, we study an interacting single-level 
quantum dot coupled to a two- level system and weakly connected to two superconducting electrodes. 
We perform a lowest-order perturbative calculation of the critical current and show that it can 
significantly change due to the two-level system. In particular, the 7r-junction behavior, generally 
present for strong interactions, can be completely suppressed. 

PACS numbers: 74.50.+r, 74.78.Na, 73.63.-b, 85.25. Cp 



I. INTRODUCTION 

The swift progress in molecular electronics achieved 
during the past decade has mostly been centered around 
a detailed understanding of charge transport through 
single- molecule junctions,—^ where quantum effects gen- 
erally turn out to be important. When two superconduct- 
ing (instead of normal) electrodes with the same chem- 
ical potential but a phase difference Lp are attached to 
the molecule, the Josephson effect implies that an equi- 
librium current I(<p) flows through the molecular junc- 
tion. Over the past decade, experiments have observed 
gate-tunable Josephson currents through nanoscale 

junctionS^^ i a»ia i n^3 1 ^,l5 1 l6d7 1 l8,l9 l 20,2l 1 22 mcmd _ 

ing out-of-equilibrium cases, and many different interest- 
ing phenomena have been uncovered. In particular, the 
current-phase relation has been measured by employing a 
superconducting quantum interference devicei 14 ' 15 ' 19 For 
weakly coupled electrodes, the current-phase relation is 3 - 



I ((p) = 7 c sin(y>), 



(1) 



with the critical current I c . 

The above questions have also been addressed 
by many theoretical works. It has been shown 
that the repulsive electron-electron (e-e) interaction 
U > 0, acting on electrons occupying the relevant 
molecular level, can have a major influence on the 
Josephson current^^^^^^^^ 2 .^ 4 * 3 ^ For 
intcrmcdiate-to-strong coupling to the electrodes, an in- 
teresting interplay between the Kondo effect and su- 
perconductivity takes place i 24 ' 28 ' 30 ' 31 ' 32 ' 35 In the present 
paper, we address the opposite limit where, for suffi- 
ciently large U, a so-called 7r-phase can be realized, with 
I c < in Eq. JT]). In the 7r-regime, tp = tt corre- 
sponds to the ground state of the system (or to a min- 
imum of the free energy for finite temperature), in con- 
trast to the usual 0-state with I c > 0, where ip = 
in the ground stated The sign change in I c arises due 
to the blocking of a direct Cooper pair exchange when 
U is large. Double occupancy on the molecular level 
is then forbidden, and the remaining allowed processes 
generate the sign change in J el 23 i 24 i 25 i 27 i 29 i 34 The most 
natural way to explain the 7r-junction behavior is by 



perturbation theory in the tunnel couplings connecting 
the molecule to the electrodes. Experimental observa- 
tions of the 7r-phase were recently reported for InAs 
nanowire dots 1 ^ and for nanotubes ; 15 ' 38 but a 7r-junction 
is also encountered in superconductor-ferromagnet- 
supcrconductor structures ! 39 ' 40 Accordingly, theoretical 
works have also analyzed spin effects in molecular mag- 
nets coupled to superconductors ^ 1 ' 42 : 43 

The impressive experimental control over supercur- 
rents through molecular junctions reviewed above implies 
that modifications of the supercurrent due to vibrational 
modes of the molecule play a significant and observable 
role i 34 ' 44 ' 45 ' 46 We have recently discussed how a two-level 
system (TLS) coupled to the dot's charge is affected by 
the Josephson current carried by Andreev states4£ For 
instance, two conformational configurations of a molecule 
may realize such a TLS degree of freedom. Experimental 
results for molecular break junctions with normal leads 
were interpreted using such models ) 48 ' 49 ' 50 ' 51 ' 52 ' 53 but the 
TLS can also be created artifically using a Coulomb- 
blockaded double dot4£ A detailed motivation for our 
model, where the Pauli matrix o~ z in TLS space couples 
to the dot's charge, and its experimental relevance has 
been given in Rcfs. H31H. While our previous work^I 
studied the Josephson-current-induced switching of the 
TLS, we here address a completely different parameter 
regime characterized by weak coupling to the electrodes, 
and focus on the Josephson current itself. We calculate 
the critical current I c in Eq. JT]) using perturbation the- 
ory in these couplings, allowing for arbitrary e-e interac- 
tion strength and TLS parameters. A similar calculation 
has been reported recently^ but for a harmonic oscilla- 
tor (phonon mode) instead of the TLS. Our predictions 
can be tested experimentally in molecular break junc- 
tions using a superconducting version of existin g 48 ' 49 ' 50 
setups. 

The remainder of this paper has the following struc- 
ture. In Sec. [TT1 we discuss the model and present the 
general perturbative result for the critical current. For 
tunnel matrix element Wq = between the two TLS 
states, the result allows for an elementary interpretation, 
which we provide in Sec. IIIII The case Wo ^ is then dis- 
cussed in Scc. lIVI followed by some conclusions in SecfVl 
Technical details related to Sec. IIIII can be found in an 
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Appendix. We mostly use units where e = h = Ub = 1- 



The Josephson current I ((f) at temperature T = (3 1 
follows from the equilibrium (imaginary-time) average, 



II. MODEL AND PERTURBATION THEORY 

We study a spin-degenerate molecular dot level with 
single-particle energy and on-site Coulomb repulsion 
U > 0, coupled to the TLS and to two standard s- 
wave BCS superconducting banks (leads). The TLS is 
characterized by the (bare) energy difference Eq of the 
two states, and by the tunnel matrix element Wo- The 
model Hamiltonian studied in this paper is motivated by 
Refs . I48ll53l where it was employed to successfully describe 
break junction experiments (with normal-state leads). It 
reads 



H = Hq + Htv 



Hi 



cads 1 



(2) 



where the coupled dot-plus-TLS part is 



Ho = — ~ ~Y° X + [ ed+ 2 CTz j + n U + Un 1 n l 

(3) 

with the occupation number n s = d\d s for dot fermion 
d s with spin s =t> 4- Note that the TLS couples with 
strength A to the dot's charge. Indeed, assuming some re- 
action coordinate X describing molecular conformations, 
the dipole coupling to the dot is oc X(n-\ + nj.), just as 
for electron-phonon couplings. 44 ' 4 ^ 48 ' 49 If the potential 
energy V(X) is bistable, the low-energy dynamics of X 
can be restricted to the lowest quantum state in each well 
and leads to Eq. © . The TLS parameters and the dipole 
coupling energy A can be defined in complete analogy to 
Refs. I48ll53l . and typical values for A in the meV range 
are expected, comparable to typical charging energies U. 
Moreover, the electron operators Cfc QS , corresponding to 
spin-s and momcntum-fc states in lead a = L/R, are 
governed by a standard BCS Hamiltonian with complex 
order parameter A L / R e zizilp ' 2 (with A L / R > 0), respec- 
tively, 



(4) 



Efe^A^cl^+h, 



where t\~ a is the (normal-state) dispersion relation. Fi- 
nally, the tunneling Hamiltonian is 



-Stun — 



where H^,J S describes tunneling of an electron with spin 
s from the dot to lead a with tunnel amplitude t\~ a , and 

r(+) - Pr(-)t 

Vas ■ 



(5) 



the reverse process is generated by Hj,J = H^ 



(6) 



where a = L/R and s =|, J, can be chosen arbitrarily by 
virtue of current conservation and spm-SU(2) invariancc, 
and T is the time-ordering operator. Equation © is 
then evaluated by lowest-order perturbation theory in 
Htun- The leading contribution is of fourth order in the 
tunnel matrix elements and can be evaluated in a similar 
manner as in Ref. H3- We assume the usual wide-band 
approximation for the leads with fe-independent tunnel 
matrix elements, and consider temperatures well below 
both BCS gaps, T < A L . R . Putting a = L and s =t, 
after some algebra, the Josephson current takes the form 
with the critical current 



T L A L dE 



T R A R dE' 



C(E,E'). 



(7) 

We define the hybridizations r Q = irpp\t a \ 2 , with 
(normal-state) density of states p R in the leads. The 
function G in Eq. ((7|) can be decomposed according to 



C{E,E')=Y J C N (E,E'), 



(8) 



N=0 



with contributions Cjv for fixed dot occupation number 
TV = nj +ni = {0,1,2}. For given TV, the two eigenener- 
gies (labeled by a = ±) of the dot-plus-TLS Hamiltonian 
H in Eq. jSJ) are 



E^ = Ne d + U5 N - 2 + 



with the scale 



The occupation probability for the state (N, a) is 



p° N = ^e-e E »(l+5 Ntl ), 



(9) 



(10) 



(11) 



where Z ensures normalization, Y^NctPn = !• With the 
propagator 



G ( (E) = 



1 



E-Z' 

we then find the contributions CV in Eq. 



(12) 
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C (E,E') 



ErP 2 r r< 7 l< 7 2&3<?4 
J 1010 



(Tl ■■■(T4 

n cr 4 /7io-iCr2CT 3 o-4 Q 



2p? T[ 



211) 



G E -2 _ B - 3 (£) G £ '2 (£') G E p _ £ -4 (-E + E') 
4_ B 'i (#) g eIS-e° 3 ( e ') G El i -El 2 (°) i 



(13) 



(71 —0*4 

+ Pi 3 G Bf3 _ SrT 4 (£') G^s _ <a (£') G E °3 _ e:i (E + E'f 



Vj_ T ^ r ^ Q ^ ^ {E) q e „. 1 _ E? {E ') G E? _ E? (E + E') 
— T^ 2 ' 73 ' 71 G E p _ E p (E) G E p _ E °i (E') G E p _ E °4 (E + E') 



(14) 



C 2 (E,E') - X ^Pl 1 T 2l2\ a3ai G El 1 - E^i^G El 1 - El 2 { E ') G El 1 - El 3 ( E + E ') 
4^2 1 1210 ^ 



(15) 



' 2 — B" 1 (-S) G B 2 °" 2 -El 3 i E ') G El 2 -E^ (°) 
I 



Here, we have used the matrix elements 

- L N 1 N 2 N 3 N 4 — ±L \ Ni N2 -ZV 3 N4 ) ' 

with the 2x2 matrices (in TLS space) 



(16) 



,± _ 1 A (So - NX)a z + W a a 
N ~ 2 I T 



iV 



For T = 0, it can be shown that Co and C2 are always 
positive, while Ci yields a negative contribution to the 
critical current. When C\ outweighs the two other terms, 
we arrive at the 7r-phase with I c < 0. 

Below, we consider identical superconductors, = 
Aji = A, and assume A > 0. It is useful to define the 
reference current scale 



In = 



TlTr 2eA 

A 2 7T 2 ft' 



(17) 



Within lowest-order perturbation theory, the hybridiza- 
tions and Tr only enter via Eq. (|17p and can thus 
be different. Equation ((TJ) provides a general but rather 
complicated expression for the critical current, even when 
considering the symmetric case A^ = A/?. In the next 
section, we will therefore first analyze the limiting case 
W = 0. 



III. NO TLS TUNNELING 

When there is no tunneling between the two TLS 
states, Wo = 0, the Hilbert space of the system can be de- 
composed into two orthogonal subspaces 7i+ © 7i_ , with 



the fixed conformational state a 
Equation ([9]) then simplifies to 

erA 
2 



± in each subspace. 



E N = 



C,J 



N + U6 N:2 



0E0 



(18) 



One thus arrives at two decoupled copies of the usual in- 
teracting dot problem (without TLS), but with a shifted 
dot level e CT = e^ + crA/2 and the "zero-point" energy shift 
—<jEo/2. As a result, the critical current I c in Eq. (|7|) 
can be written as a weighted sum of the partial critical 
currents I c (£<r) through an interacting dot level (without 
TLS) at energy e CT , 



<y=± 



(19) 



where p a — ^2 N p" N with Eqs. (HB) and pj denotes the 
probability for realizing the conformational state a. The 
current I c (e) has already been calculated in Ref. [H (in 
the absence of phonons), and has been reproduced here. 
In order to keep the paper self-contained, we explicitly 
specify it in the Appendix. 

In order to establish the relevant energy scales deter- 
mining the phase diagram, we now take the T = limit. 
Then the probabilities (jlip simplify to p a N = Sj^j^dug-, 
where E°^ = mincer) (E^) is the ground-state energy of 
Hq for Wo = 0. Depending on the system parameters, 
the ground state then realizes the dot occupation num- 
ber N and the TLS state a. The different regions (N,a) 
in the Eo — td plane are shown in the phase diagram in 
Fig. Q] The corresponding critical current in each of these 
regions is then simply given by I c = I c (e s ). 
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FIG. 1: Ground-state phase diagram in the Eo — €d plane 
for Wo = 0. Different regions (TV, a) are labelled according 
to the ground-state dot occupation number TV = 0, 1,2 and 
the conformational state a = ±. Dark areas correspond to 
7r-junction behavior. The charge-degeneracy line = —U/2 
is indicated as dashed line. Main panel: A < U. Inset: A > U, 
where no it -junction behavior is possible for U < Eo < 2A — U. 

By analyzing the dependence of the ground-state en- 
ergy on the system parameters, one can always (even for 
Wo ^ 0) write the function C(E,E') in Eq. J7J as 

C(E,E') = Q(^-e d )C 2 (20) 
+ Q(e d - £_) 0(£+ - e d ) C x + 6(e d - £+) C , 

where is the Hcaviside function and the energies £± = 
£±(U, A, E ) are the boundaries enclosing the 7r-phase re- 
gion with N = 1, i.e., (£_) denotes the boundary 
between the N ~ and N = 1 (the TV = 1 and TV = 2) 
regions, see Fig. [1] Explicit results for £± follow from 
Eq. Ull) for VF = 0. For E < (£b > 2A) and arbitrary 
TV, the ground state is realized when a = — (a = +), 
leading to £ + = A/2 (£ + = —A/2). In both cases, the 
other boundary energy follows as £ £ . U. In the 
intermediate cases, 

with£ = \{X-U-E ), wc find for 

< E < A, 

£+ = max(A/2 - Sq, = min(A/2 - U, fo), (21) 

while for A < So < 2A, we obtain 

£+ = max(-A/2, £ ), = min(C , A/2 + 2f ). (22) 

These results for t;± are summarized in Fig. [1] Remark- 
ably, in the Eq — plane, the phase diagram is inversion- 
symmetric with respect to the point (Eq = A, ta = 
— U/2). Furthermore, we observe that for many choices 
of Eo, one can switch the TLS between the a = ± states 
by varying e^, sec Fig. [T] 




o- 



-2-10 1 

FIG. 2: (Color online) Ground-state critical current I c as a 
function of for Wo — 0. / c is given in units of 7o, see 
Eq. (|17|l . In all figures, the energy scale is set by A = 1. 
Dashed (red), dotted (blue) and solid (black) curves represent 
the partial critical currents 7 c (e+), 7 c (e_), and the realized 
critical current I c , respectively. Main panel: Eo = 0.8, A = 
0.6, U = 1, such that > This corresponds to the 7r- 
phase region with A < Eo < 2A in the main panel of Fig. [T] 
Inset: Eo = 0.6, A = 0.8, U = 0.5, where < and no tt- 
junction behavior is possible. This corresponds to U < Eo < 
2A — U, see inset of Fig. [T] 



immediately conclude that the 7r-junction regime (where 
N = 1) can exist only when £ + > This condition is 
always met away from the window < Eq < 2 A. How- 
ever, inside that window, Eqs. (f2Tj) and |22)) imply that 
for sufficiently strong dot-TLS coupling, A > U, the ir- 
phasc may disappear completely. Indeed, for U < Eo < 
2 A — C7, no 7r-phasc is possible for any value of e d once 
A exceeds U . The resulting ground-state critical current 
is shown as a function of the dot level e d for two typical 
parameter sets in Fig. [51 The inset shows a case where 
the 7r-phase has been removed by a strong coupling of the 
interacting dot to the TLS. The above discussion shows 
that the 7r-junction regime is very sensitive to the pres- 
ence of a strongly coupled TLS. 



IV. FINITE TLS TUNNELING 

Next we address the case of finite TLS tunneling, Wo ^ 
0. Due to the a x term in Ho, the critical current cannot 
be written anymore as a weighted sum, see Eq. p^j) . and 
no abrupt switching of the TLS happens when changing 
the system parameters. Nevertheless, we now show that 
the size and even the existence of the 7r-phase region still 
sensitively depend on the TLS coupling strength (and on 
the other system parameters). In particular, the 7r-phasc 
can again be completely suppressed for strong A. 



We now notice that Eq. (|20[) implies the same decom- 
position for the critical current ([7]). Wc can therefore 
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FIG. 3: (Color online) Phase diagram and boundary energies 
£± enclosing the 7r-phase for finite Wo- Main figure: A = 0.4 
and U = 0.5, where a 7r-phase is present; Wo = 0,0.2 and 5, 
for solid (blue), dashed (black) and dash-dotted (red) curves, 
respectively. Inset: A = 0.7 and U = 0.5, where the 7r-phase 
vanishes; Wo = 0,0.3 and 3, for solid (blue), dashed (black) 
and dash-dotted (red) curves, respectively. 

For finite Wo, the ground-state critical current is 
obtained from Eq. (|2"0|) , where the Cm are given by 
Eqs. l[T3 j) -([TB ]l and the 7r-phase border energies £± are 
replaced by 

f+ = 5(*i-*o), £- = i(* a -*i-2l0. (23) 

The 3>tv are defined in Eq. (fT0|) . Compared to the Wo = 
case in Fig. [1] the phase diagram boundaries now have a 
smooth (smeared) shape due to the TLS tunneling. Nev- 
ertheless, the critical current changes sign abruptly when 
the system parameters are tuned across such a boundary. 
The energies (|23|) are shown in Fig. [3] for various values 
of Wq in the Eq — td plane. In between the and in- 
curves, the 7r-phase is realized. From the inset of Fig. O 
we indeed confirm that the 7r-phase can again be absent 
within a suitable parameter window. Just as for Wo = 0, 
the 7r-phase vanishes for £ + < £_ , and the transition be- 
tween left and right 0-phase occurs at £ = (£+ + £-)/2. 
For \Eo\ 3> max(A, Wo), we effectively recover the phase 
diagram for Wo = 0, since the TLS predominantly occu- 
pies a fixed conformational state. 




FIG. 4: (Color online) Listline plots of the Wo / ground- 
state critical current I c (in units of Io) in the Eo — £d plane, 
with A = 0.6 and U = 0.5. The boundaries £± enclosing the 
7r-phase, see also Fig. [3] are indicated as solid (blue) curves. 
Top panel: Small tunnel amplitude, Wo = 0.2. Bottom panel: 
Large tunnel amplitude, Wo = 3. 



realized in a side-coupled double-dot system^ 7 - Finally 
we note that, unlike for Wo = 0, the perturbative result 
for the critical current diverges at the point where the 
7r-phase vanishes, i.e., for = £. This divergence is 
an artefact of perturbation theory and is caused by the 
appearance of the factor G E - _ E - (0) = (e<j — £) _1 in 
Eqs. P]) and ([I5]). 



V. CONCLUSIONS 



The corresponding critical current I c is shown in Fig. [4] 
for both a small and a very large TLS tunnel matrix 
element Wo. In the limit of large Wo 3> max(A, |-Eo|)j see 
lower panel in Fig. 0] the dot and the TLS are effectively 
decoupled, since (<r z ) ~ and (cr x ) ~ sgn(Wo). While 
this limit is unrealistic for molecular junctions, it may be 



In this paper, we have presented a perturbative cal- 
culation of the critical Josephson current, I Cl through an 
interacting single- level molecular junction side-coupled to 
a two-level system (TLS). Such a TLS is a simple model 
for a bistable conformational degree of freedom, and has 
previously been introduced in the literature! 47 ! 48 ' 53 Our 
perturbative calculation assumes very weak coupling to 
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attached superconducting reservoirs. The ground-state 
critical current can then be computed exactly for other- 
wise arbitrary parameters. Our main finding is that the 
7r-phase with I c < is quite sensitive to the presence of 
the TLS. In particular, for strong coupling A of the molec- 
ular level to the TLS as compared to the Coulomb energy 
U on the level, the 7r-phase can disappear altogether. 
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where 



C%(E,E') = ir N c° N (E,E') 



Pn = 



Z, 



L e -/^ ( i + (W)) 



with Z a such that ^2 N Ppj 
given by 



= 1. Moreover, the c a N are 
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c%(E,E') = 



1 



(E + e a )(E> + e a ) 



1 



E + E' 2e a + U 



APPENDIX A: PARTIAL CRITICAL CURRENTS 

In this Appendix, we provide the partial critical cur- 
rent ic(e CT ) which appears in the calculation for W = 0, 
see Sec. IIIII In the absence of TLS tunneling, the matrix 
elements (fl6j) simplify to 



rpO 1(72^3^4 

1 N 1 N 2 N 3 N i 



cl (E,E') 



1 



E + E' 



1 



{E-e a ){E + e a + U) 
1/2 



(E> -e a )(E> + e a + U) (E - e a ){E> - e a ) 
1/2 



{E + e a + U)(E' + e a + U) 



where cr, = er.;Sgn(A^A — Eq). We now rename a — > a to 
denote the conformational state (eigenstate of a z ). 

The partial current / c (e<r) corresponding to fixed con- 
formational state a = ± is then given by 



Ic{e a ) 
h 



^ 3 E 

N 



dEdE'C^(E,E') 
A y/(E 2 - A 2 ) (£' 2 - A 2 



4{E,E') 



1 



{E - e a - U){E' - e a - U) 
1 2 



E + E> 2e a + U 
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